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RELIABILITY ANALYSIS OL SYSTEMS 
SUBJECT TO FIRST-PASSAGE FAILURE 


INTRODUCTION 

An obvious goal of reliability analysis is the avoidance of system failure. However, 
it is generally recognized that it is often not feasible to design a practical or useful system 
for which failure is impossible. Thus it is necessary to use techniques that estimate the 
likelihood of failure based on modeling the uncertainty about such items as the demands 
on and capacities of various elements in the system. This usually involves the use of 
probability theory, and a design is considered acceptable if it has a sufficiently small 
probability of failure. In some applications, the analysis uses some rule-of-thumb, such as 
three-sigma, six-sigma, etc., rather than numerical estimates of the probability of failure, 
but these approaches are based on known correlation between the rule-of-thumb and 
estimated probability of failure in a random-variable model. 

A random variable is a mathematical concept to model the uncertainty about any 
quantity that has more than one possible numerical value. That is, it models the likelihood 
that the observed outcome will be any particular one of the possible outcomes. The 
situation becomes more complicated when the time of likely failure is also an uncertain 
quantity, such as when the demand on the system varies randomly with the passage of 
time. Then, one must use a method that models not only the uncertainty about quantities 
at any particular time, but also the uncertainty about how these quantities vary from one 
time to the next. The probabilistic model that combines these two types of uncertainty is 
called a stochastic process (or random process). Thus, the uncertain stress at one instant 
of time t will be classified as a random variable X(t ) but the uncertain history of stress 
over a range of time values will be a stochastic process {X(t)} . The practice of denoting 
a stochastic process by putting the notation for the associated random variables in braces 
will be used to indicate that the stochastic process is a family of random variables — one 
for each t value. 

Study of first-passage failure is appropriate when the system of interest is 
considered to have performed unsatisfactorily if some measure of response has ever 
reached some critical value. Thus, for example, first-passage failure might be considered 
to have occurred if the stress or strain at some critical location has ever exceeded the 
yield level or if a particular displacement has exceeded a value resulting in collision with 
some other element of the system, possibly causing damage. 

To calculate the probability of first-passage failure during a given time interval 
[0,f], one needs tenns such as P[X(t)> x crilical \ in which X(t) is defined as the 
maximum of X(s) over the set of values 0 < s < t . For a given value of t, the first-passage 
reliability can be formulated in terms of the random variable X(t), but finding (or 
estimating) the probability distribution of X(t) requires a significant amount of 
information about the stochastic process \X(s ) : 0 < s < t) . 



A fundamental difference between a random variable and a stochastic process 
relates to what constitutes one observation of an outcome. An outcome for a scalar 
random variable is a scalar. An outcome for a vector random variable is a vector. But an 
outcome for a scalar stochastic process is a time history. Furthennore, most stochastic 
processes have the ergodic property, which assures that any particular time history x{t) 
must eventually approach every possible value of the stochastic process {X(t) \ . While it 
is relatively rare that an observation of a random variable will exceed the mean plus four 
standard deviations of the variable, for example, it is expected that any one observation of 
the associated process will eventually reach that four-sigma level if the observation time 
is sufficiently long. Thus, the length of time f is a crucial variable in first-passage 
analysis. 

It may be noted that there is a sort of “orthogonality” between the idea that a 
stochastic process is characterized by an “ensemble” containing all possible time histories 
and the idea that the process is a family of random variables. A time history is a single 
observation including many values of t. Many observations at a single value of t give a 
statistical sample of the random variable X(t). In most practical problems, it is not 
feasible to describe a stochastic process in terms of the probability of occurrence of 
particular time histories, but it is always possible to characterize the process by using 
information about the joint probability distribution of the random variables of which it is 
composed. 


PROPERTIES OF STOCHASTIC PROCESSES 

Almost every property of a stochastic process can be viewed as a generalization of a 
corresponding property of a random variable. Of course, all the information about the 
probabilities of a single random variable X(t) can be encapsulated in a cumulative 
distribution function 


F x(t)( u ) = p i x (t) ^ u\ 


( 1 ) 


or by its derivative, which is the probability density function 


, , dF x(l) (u) 


.. P[u < X(t) < u + Au\ 

Inn — - 

A u — ^ 0 /l II 


( 2 ) 


Note that both the functions F x ^(u) and p X (t )( u ) must be defined on the set of all 

possible values of X{t). In order to generalize (1) and (2) to describe a stochastic process 
{X(t)} one must consider the probabilities of events that are the intersections of 
occurrences at multiple values of t. In particular, a complete description requires 
knowledge of 


F X(t l ) - X(t„)( u \’'"-> u n) = F \- X ( t \)- u W’X(tn)- u tr\ 


( 3 ) 


2 



or 


,_d lF X{t x )-X{t n )(. u b"-’ U n) 
Px( h )• ■ -X(t„ )K u W’ U n) = 


du v --du n ^ 

P[ui < X(t\) < u + Au l ,---,u n <X(t n )< u + Au n \ 

Au [ — >0,“\Au n — >0 Au^''‘Al1ji 


and one of these functions must be known for all possible choices of n, t x ,---,t n , and 

In many situations it is possible, and much easier, to use an analysis that requires 
less infonnation than is given by (l)-(4). In particular, the idea of the mean and variance 
values of random variable are generalized to a mean-value function 


Mxtt) = E[X(t)\ = J x u p xu } (u) du 


( 5 ) 


and a covariance function 

^ xx (P s ) = E([X(t) - Px(t)][X(s) - 

r oo r co 

= L x ~Mx(t) ][«2 -Mxis)} Px{t)X{s)( u b u 2) du \ du 2 


( 6 ) 


Note, in particular, that the covariance function is defined on a two-dimensional set 
of (7,5) values. It includes the much simpler variance of all the random variables, since 

a X(t) = Kxx(td), but it also includes crucial information about the correlation of the two 
random variables X(t) and X(s) . In some situations it is also necessary to consider 
covariance functions involving two different stochastic processes. One such example is 
the cross-covariance of a process and its time-derivative 


K xx (t,s) = E([X(t) - fi x (0 IX (s) -tu x (s)])= (7) 

os 

The property of stationarity (or homogeneity) of a stochastic process { X(t)\ always 
refers to some aspect of the description of the process being unchanged by any arbitrary 
shift along the t axis. There are many types of stationarity depending on which 
characteristic of the process has this property of being invariant under a time shift. 

The simplest type of stationarity involves only invariance of the mean value 
function for the process. In particular, {X(t)\ is mean-value stationary if 
p x {t + r) = ju x (t ) for any value of the time-shift parameter r. Clearly this can be true only 
if ju x (t) is the same for all t values, so one can say that { X(t)\ is mean-value stationary 
if ju x (t) = ju x , in which the absence of a t argument on the right-hand side conveys the 
information that the mean value is independent of time. Although the notation /j x is the 
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same as for the mean value of a random variable, it is here used to refer to the mean value 
function of a stochastic process. Of course, having the mean value be independent of t 
does not imply that the {X(t)\ random variables are all the same at different values of t 
or that the probability distributions of these random variables are all the same — only that 
they all have the same mean value. 

Similarly, one can define a more rigorous stationarity related to the covariance 
function. The process {X(t)} is said to be covariant stationary if the covariance function 
is invariant under a time shift so that K-^it + r,s + r)=K xx (t,s) for any value of r. 
Choosing r = -s, then gives K xx (t,s) = K xx (t-s, 0), showing that the stationary 
covariance function depends on only one time argument — the difference between the two 
times of interest. Thus, one can define a covariance function of only one time argument 
as G jQf(z') = K z - , t) so that 

Kxx(t,s) = Gxx(t - s ) for any values of t and 5 (8) 

Similarly, two stochastic processes {X(t)\ and \Y(t)\ are jointly covariant stationary if 
K xy'( l + r - s + r ) = X x Y (Ys), and it is convenient to define G X7 (r) = z,t) so that 

Kxy(t,s ) = Gxy(t - s ) for any values of t and 5 (9) 

Although a given stochastic process may simultaneously have various types of 
moment stationarity, this is not necessary. In particular, a process may be mean-value 
stationary and covariant nonstationary, and another process may be mean-value 
nonstationary and covariant stationary. It is common, however, to have situations in 
which both the mean-value and the covariance are stationary. 

There are also forms of stationarity that are not defined in terms of moment 
functions. Rather, they are defined in terms of probability distributions being invariant 
under a time shift. The general relationship is that (X(t)} is /th-order stationary if 

PX(t ] +r)--X(tj+r ) (u h ---,Uj) = Px (h )...X(t j )(u h -",Uj) ( 10 ) 

for all values of {ti,---,tj,u h ---,Uj} and the shift parameter r. This includes, as special 
cases, (X(t)} being first-order stationary if Px(t+r)( u ) = Px(t)( u ) anc * second-order 
stationary if p X (t 1+ r)X(t 2 +r)( u b u 2) = Px( tl )X(t 2 )( u b u 2) for all values of {t h t 2 ,u h u 2 } and 

the shift parameter r. First-order stationarity always implies mean-value stationarity, and 
second-order stationarity implies both mean-value and covariance stationarity, In some 
situations (such as the well-known Gaussian process) mean-value and covariant 
stationarity also imply second-order stationarity, but this is not true in general. 

The most restrictive type of stationarity is called strict stationarity. The process 
(X(t)} is said to be strictly stationary if it is /th-order stationary for any value of j. This 
implies that any order probability density function has time shift invariance and any order 
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moment function has time shift invariance. In common usage, one can also say that the 
word stationary without qualifier simply means that all moments and/or probability 
distributions being used in the given problem are invariant under a time shift. 


An ergodic stochastic process is one for which any infinitely long time history is 
representative of the entire process. For example, any expected value of some function of 
the process can be written as the limit of a time-average obtained from one time history 
of the process. Thus, 

Mx = E [ x (t)\ = lim ^ \ T x(t)dt (11) 

T — »co 7 u 


and 

G xx (t) = E([ x (t+ T)-jU x ][X(t)-jU x ]) 

If T 

= lim T Jn W> + T )-Mx][x(t)-Mx]dt 
T — »C0 1 U 


( 12 ) 


with x(t ) representing any one time history of { X(t)\ in both (11) and (12). Note 
that the process is required to be mean-value and covariant stationary in (11) and (12), 
respectively, since the infinite time average could not possibly converge to a function of 
t. Similarly, strict ergodicity implies that the probability distribution can also be found 
from one time history: 

F XU)( U )= E ( u [ u ~ x ( t )])= 11 ™^-\ TI ~U[u-x(t)\dt (13) 

7— »oo L u 


in which the unit step function U(-) is defined to be zero when its argument is 
negative and unity otherwise. Thus, the integral gives exactly the amount of time for 
which x(t) < u within the interval [0,7]. This condition that the probability that X(t) < u 
at any time t is the same as the fraction of the time that any time history x(t) is less than 
or equal to u is the basis of the statement in the introduction that any time history of an 
ergodic process must eventually approach every one of its possible values. Although it is 
not necessary for all processes to be ergodic, it is common practice to assume ergodicity 
unless there is some obvious physical reason why that would be inappropriate for the 
problem of interest. 

In the remainder of this presentation, the process |X(t)} will be considered to be 
mean-zero. The primary motivation for doing this is that it simplifies the presentation 
without losing any generality in almost all instances. In order to analyze a process {F(?)} 
that is not mean-zero, one simply defines X(t) = Y(t)~ fj. Y (t). Furthennore, this 
transfonnation is generally not difficult, since the evaluation or estimation of the mean- 
value function is generally much easier than evaluation of other quantities such as 
covariance. This is illustrated in (11) and (12) for the stationary ergodic situation, and it 
is also generally true in the analysis of the response of dynamical systems [e.g., Lutes and 
Sarkani, 2004], 
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The Gaussian stochastic process occupies a very special place in analysis of 
practical problems. It has been found that it both gives a reasonable model of many 
physical processes and has some very desirable mathematical characteristics. The 
definition of a Gaussian stochastic process {X(t)} is simply that any finite set 
{X(t l ),X(t 2 ),- ■ • ,X(t n )} of random variables from that process has the jointly Gaussian 
distribution. The general form for the joint probability density function for a vector V of 
jointly Gaussian mean-zero components is 


/>v( u ) = - 


1 


{2n) nl2 K 


il/2 


exp 


w 


■iu r K-‘ 
2 


w 


u 


(14) 


T 

which is a function of the square covariance matrix K vv = E(X V). Applying this 
model to the Gaussian process, then, only involves using the vector 
X = [X( tl ),X(t 2 ),---,X(t n )] T . 

The jointly Gaussian distribution, of course, implies that each of the random 
components is Gaussian (nonnal), but it also specifies much more about the 
interrelationships among the components. One of the very desirable features of the mean- 
zero Gaussian process is that the distribution depends only on the covariance function 
Kxx{t,s) or Gxx{z), which provides all components of Kyy for the process. Among the 
unusual features of the distribution are the facts that a mean-zero covariant stationary 
Gaussian process is always strictly stationary, and that a mean-zero covariant ergodic 
Gaussian process is always strictly ergodic. 

Particularly important for many applications is the fact that any linear combination 
of jointly Gaussian random variables is itself Gaussian, and it is also jointly Gaussian 
with other such linear combinations. This ensures that linear operations on a Gaussian 
stochastic process always yield other Gaussian processes and that these new processes are 
jointly Gaussian with the original process. Such linear operations include input-output 
relationships for any linear dynamical system and calculus operations such as finding the 
derivative process { X(t)} corresponding to a given process {X(t)\ or integrating to find 
the {X(t)} corresponding to a given \X(t )) . It should also be noted, though, that there 
are some processes for which the random variables X{t) are all individually Gaussian, 
but not jointly Gaussian. These are not Gaussian processes, and applying a linear- 
combination operation to them generally does not yield a Gaussian output. 


FREQUENCY, BANDWIDTH, AND AMPLITUDE 

The Fourier transfonn provides a classical method for decomposing a time history 
into its frequency components, and this is often very useful in both computing and 
interpreting the results of practical models. When applied to a stochastic process (X(t)} 
the Fourier transform introduces a new stochastic process {X(co)} defined by 
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( 15 ) 


X(aj) = — j°° X{t)e~ icot dt 
In ~°° 

The original stochastic process can then be recaptured by using the inverse 
transfonn 

X(t)= J ^ X(co)e t0H doj (16) 

Note that { X(co)} gives a complex random variable X(co) for every possible value 
of co, just as {X(t)} gives a real random variable X(t ) for every possible value of t. 

From (15) one can write E[X{co l )X{co 2 )\ as a double Fourier transform of G xx { r). 
A major difficulty in applying the Fourier transform procedure to many problems of 
interest, though, is the fact that E[X(( 0 \)X(a> 2 )\ does not exist for all values of frequency 
when {X(t)} is a stationary stochastic process. That is, the double integral to evaluate 
E[X((Q\)X(<x> 2 )\ ma y be infinite. It is important, though, to have a usable form of 
E[X{co l )X{co 2 )\- To this end, the Fourier transfonn procedure has been modified as 
follows. The first step is to truncate the Fourier transfonn integral in (15) to be over a 
time interval of finite length T, which is found to give an approximation of 
E[X(cq\)X(<x> 2 )\ that grows in proportion to T. This is then normalized to give a finite 
limit, and is defined as the spectral density function of the process: 

y *,][T d, 

= fZ'ia bm i", 2 2 G «('l - '' t " 2 ' 2) <*1 d ‘2 

For any stationary process, though, it has been shown that (17) gives the spectral 
density as exactly the Fourier transform of the covariance function for the original 
process: 

Sxx{®)=-^\_ co G X x{i:)e lC0T dr (18) 

Furthermore, this implies that the inverse must also be true: 

Gxr(r) = \Zjxx(oS)e l0} r dco (19) 

The original idea that E[X{co{)X{( 02 )\ is a double Fourier transform of G xx ( z) 

gives 

G xx (z) = j' ; J ^ EiXUo, )X(co 2 )]e l(a> ' +( '^- 1 dco x dco 2 (20) 
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and this agrees with the expression in (19) if E[X(a) l )X(co 2 )] = Sxx(<^> + <v 2 )’ i n 
which <5(-) denotes the Dirac delta function. This degenerate fonn has been found to be 
useful in some calculations. 

The spectral density function defined in (17), and commonly computed from (18), 
has many desirable properties. Most importantly, it does give a complete frequency 
decomposition of the covariance properties of a stationary {X(t)} . In addition, it is 
always a real, nonnegative even function of co. Furthermore, setting z equal to zero in 
(19) shows that the variance <j x = G xx (0) is simply the area under the S xx (co) curve. It 
should be noted that there are other forms of spectral density that differ from (18) only by 
a multiplying factor. One such variation is found by using a cosine Fourier transform in 
place of the exponential transform used here, and obtaining a spectral density that is 
twice as large as that of (1 7)-( 19), but exists only for co > 0. A less frequent variation is to 
use frequency measured in Hz (cycles per second) rather than radians per second in 
saying that cr x is the area under the spectral density curve, which requires a modification 
of spectral density by a factor of 2 n. 

The cross-spectral density between two processes, such as |X(t)} and |X(t)} is 
defined in a parallel way, and (18) and (19) generalize to give the cross-spectral density 
and the cross-covariance function as a Fourier transform pair. Taking derivatives of (18) 
and (19) yields simple results for the cross-spectral density of {X(t)} and \X(t)\ , the 
spectral density of |X(t)}, etc. In particular, S xx (co) = -icoS xx (co) and 

S xx ( 0J ) = (0 ^xx^ (0 ) ■ The results for higher order derivative processes can be written as 
S v ( j)v (k){o)) = (~\) k (i) J+k ( 0 J+k S xx (co), in which \X {j) (t) } represents the /th-order 

derivative with respect to t. This includes the special case of S x(j)x(J) (co) = co 2j S xx (co) 
for the spectral density function of the /th derivative process. 

A process {X(t)} is said to be narrowband if the spectral density S xx (co) is very 
small except within a narrow band of frequencies. Because S xx (co) is an even function, 
this really means that the band of significant frequencies appears both for positive and 
negative co values. One way to state this narrowband condition is to say that S xx {co) « 0 
unless \co\~ co c for some given characteristic frequency co c . The limiting narrowband 
process has the spectral density S^ico) = (cj x I2)[8(a>+ a> c ) + S(co- co c )]. The 

'y 

corresponding covariance function is G xx (t) = <j x cos(o> c t). Furthennore, the time 
dependence of the process is of the form X{t) = A eos[<» c t + 6\ with the amplitude and 
phase A and 9, respectively, being random variables. That is, any time history is a pure 
cosine wave with fixed amplitude and frequency. This situation, of course, is not often 
encountered in practice, but the harmonic nature of X{t) and G xx ( z) is approximated by 
all narrowband processes. The time histories of X(t) for a narrowband process generally 
are of the form of cosine waves with slowly varying amplitude and phase, as shown in 
Figure 1. Similarly, the corresponding G xx ( z) has the form of a cosine wave with 
amplitude (or envelope) that decays slowly with increasing | r |. 




At the opposite extreme from a narrowband process is the white noise (or delta- 
correlated) process for which all Fourier components contribute equally, such that the 
spectral density is the same for all co values. For a process {X(t)} of this type with 
S xx ((o) = Sq, the covariance function is the degenerate case of Gfyr(r) = 2/rS 0 8{z). The 
time histories of any delta-correlated process, of course, must be extremely erratic, since 
they are uncorrelated at times separated by any finite value r. Any practical broadband 
process has a spectral density that varies slowly with co, a covariance function that tends 
rapidly to zero with increasing r |, and very erratic time histories. 

Many stochastic processes encountered in practice are not near either the 
narrowband or the broadband limits. Some of these may be considered to be between 
narrowband and broadband, while others have distinctly different forms, such as a 
multimodal spectral density containing a sequence of narrow peaks at distinct values of 
co. In this case, many but not all frequencies make significant contributions to the time 
histories. One physical phenomenon that leads to narrowband or multimodal spectral 
densities relates to the deterministic idea of resonance in dynamical systems. If a system 
has a resonant frequency, then exciting that system with a broadband process will lead to 
a response for which the spectral density has a high peak at the resonant frequency. In the 
same way, multiple resonant frequencies in the system can lead to multimodal spectral 
densities. In essence, any dynamical system can be regarded as a filter. Some frequencies 
components of the excitation are attenuated because the system is not sensitive to those 
frequencies, while other components are amplified due to resonance-like behavior. For a 
linear system or filter with stochastic input | F(/)| and output {A(/)j there exists a 
complex transfer function II XY ((o) such that the Fourier transforms satisfy 

X(co) = Hxy(o))Y(co). This then gives S x (co) =|// xy (<u) |~ S Y (oo) as the description of 
the filtering effect on the spectral density. 

The concepts of characteristic frequency, amplitude and phase that were introduced 
in discussing narrowband processes can be usefully generalized to provide an alternate 
characterization of any mean-zero stationary stochastic process. In particular, let 

X{t) = A{t)cos[co c t+ Oit)] (21) 
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The definition of co c should ideally be chosen such that E[6(t)\ = 0 , because a 
nonzero value for E[6(t )] should, in itself, be a contribution to the characteristic 
frequency of the \X(t)\ process. Next the amplitude and phase are defined as 

A(t) = [X 2 (t) + Z 2 (t)] m (22) 

and 

0(t) = -tan _1 \Z(t)/X(t)\-co c t (23) 

in which { Z(t)) is a process for which every frequency component is 90° out of phase 
with the corresponding component of |X(f)} . Specifically, this is done by choosing 
Z(co) = ig(co)X(co) for some odd real function g(-) with g(co c ) = 1. This definition of the 
amplitude and phase is consistent with the intuitive ideas introduced for a narrowband 
process for which X(co) is nearly zero except in the neighborhood of ±co c , and there are 
many g(-) functions satisfying the necessary conditions. One of the popular choices is 
g(co ) = sgn(<z>), which gives Z(t) as the Hilbert transfonn of X(t). The amplitude defined 
in this way was used by Cramer and Leadbetter [1967] and is commonly associated with 
their names. It will be denoted as A x (t ) . Another commonly used amplitude follows from 

using g(co) = co I co c , which gives Z(t ) = X(t)i co c . This amplitude, which will be denoted 
as A 2 (t), is often referred to as the energy-based amplitude, since it can be related to the 
sum of potential and kinetic energy in some simple oscillators. The co c values 
corresponding to the Cramer/Leadbetter and energy-based amplitude are 2 , / A () and 

i / o 

(Z 2 /Z 0 ) , respectively, in which the spectral moment is defined as 

f °0 |* 00 

Zj= J | coy Sxxt 00 ) dco — 2J q co S xx (co)dco (2d-) 

The term bandwidth parameter is used to refer to any characteristic of a stochastic 

process that gives an indication of how strongly the process resembles the narrowband or 
broadband limit. There are many such parameters, but only two will be mentioned here. 
One natural way to introduce this idea is to consider the rates of change of A(t) or 6(t), 
as defined in (22) and (23). In particular, these quantities should change slowly for a 
narrowband process, but more rapidly for a broadband process. For the Cramer/ 
Leadbetter amplitude it can be shown [e.g., Lutes and Sarkani, 2004] that the rates of 
change of A x (t) and of the corresponding phase are governed by the parameter 

a x =Z x !(Z () A 2 ) m (25) 

Similar results apply for the energy-based amplitude and phase, but the controlling 
parameter is 

a 2 = /L 2 /(/L 0 /L 4 ) 1/2 (26) 
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It may be noted that the even spectral moments have direct physical meaning since 
<7 x ~ 4> > a x = ^2 > a x = ^4 ’ etc - This shows that the value of a 2 can be found directly 

from these variance quantities, which is something of an advantage over a 1 . On the other 
hand, a 2 depends rather heavily on the high-frequency tail of the spectral density, which 
is sometimes a disadvantage in comparison to a 1 . Both a 1 and a 2 are defined such that 
they are in the interval [0, 1] and tend to unity for narrowband processes. 

Note that if {X(t)} is a Gaussian process, then {X(a>)} , {Z(<x>)},and {Z(t)} are all 
jointly Gaussian processes, because they are obtained from linear operations on |X(t)} . 
The Gaussian property then allows evaluation of additional properties of the amplitude 
and phase. It can be shown that both the Cramer/Leadbetter and energy-based definitions 
give the random variable Z(t) to be independent of X(t) for the same t value and to have 
the same variance. This then leads to A(t ) having the common Rayleigh distribution: 

F a (u) = \- e - “ 2 /(2c7 ^ } p A (u) = {u/ax)e~ u2l(2a ^ ) (27) 

The corresponding phase angle 6(t) is independent of A(t) for the same t value and 
is unifonnly distributed on [0,2;r]. These simple distributions are sometimes quite useful 
in evaluating properties of (X(t)} from the amplitude and phase form of (21). 


CROSSING PROBABILITY AND CROSSING RATE 

The reliability of a system is closely related to the concept of level-crossings. This 
is particularly true for first-passage failure, in which the system is considered to fail only 
when a particular stress process or displacement (X(f)} reaches a critical level b. In 
particular, failure cannot then occur before X(t) crosses b. Similarly, a crossing of the 
level zero by the derivative X(t) coincides with a local maximum or minimum of X(t ) , 
and failure will not occur when all these extrema are below b. The following paragraphs 
review important results related to these crossing probabilities, with special emphasis on 
the stationary Gaussian process so commonly used in practical applications. Rigorous 
analytical results will be presented for the expected rate of occurrence of such crossings 
for this process, and for the corresponding probability distribution of local maxima. 
These results, however, cannot provide an analytical solution for the first-passage 
probability so further approximations will be introduced for that purpose. 

An upcrossing at time t of the level it by the process \X(t)\ is simply the 
intersection of the events X(t) = u and X(t) > 0 . Similarly, a downcrossing occurs 
whenever X(t) = u with X(t) < 0 . These ideas are illustrated in Figure 2. The 
probabilities of such crossings during a small time increment will be written as 

Vx(u,t)z 1 1 = /’(upcrossing in \t,t + A t]) (28) 
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and the corresponding relationship between a downcrossing and a term v~ x (u,t). Thus, 
v x (u,t) and v x (u,t) are the limits of /^crossing in [?, t+At])/At as At tends to zero. 
Note that v x (u,t) and v x (u,t ) have units of crossings per unit time, and are exactly the 
expected rate of crossings of the level X=u. The terms upcrossing rate and 
downcrossing rate will be used to refer to these quantities, even though the interest here is 
more focused on the probability of occurrence of the crossing than on the rate. 

It is easily seen that there can only be an upcrossing of the level u within the 
interval [/, t+At] if X(t) at the beginning of the interval is less than u, but close to u, 
and has a positive derivative. Inasmuch as A t is infinitesimal, one may consider the 
derivative to be constant at the value X(t) throughout the time interval and conclude that 
there will be an upcrossing within the interval if u - X(t) At < X(t) <u . This event is 

shown shaded on the space of possible values of X( t) and X (?) in Figure 3. The 
probability of this event can now be found by integrating the joint probability density of 
X(t) and X(t) over the shaded region. However, the fact that At is infinitesimal can be 
used again to argue that P X(t ) X ^ t ){w,v) « P X (t)X(t)^ u,v ) for u-v At <w <u, giving: 

poo poo 

y + x(p,t) = J o vp X(t)m (u,v)dv = p X( 0 (m)J o vp m [v I X(t) = u]dv (29) 

The second form in (29) has been obtained simply by writing, the joint probability 
density function as the product of a marginal and a conditional density function. 
Similarly, the downcrossing rate is 

v x {ii,t) = -\ Q ^ v P X {t)X{t)^ u ^ dv = Px(t)( u ) J_j v \Px{f)W(f) = u ¥ v (3°) 

The results in (29) and (30) are rigorous for any stochastic process for which the 
integrals exist. It should be noted, though, that the values obtained might be infinite if the 
conditional probability density function of X{t) does not decay sufficiently rapidly. For 

example, if E[X{t)\X{t) = u] does not exist then v x (u ,t) may be infinite. For a 
stationary { X(t)} process these rates, of course, are independent of time, and will be 
denoted simply as v x (u) and v x (u) . 
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X(t) 


X(t) = u-X(t 



X(t) = u 


X(t) 


Figure 3. Phase Diagram showing the Event of Upcrossing 


For any covariant stationary process, the random variables X(t ) and X(t) are 


uncorrelated. In the special case of a Gaussian process, they are then uncorrelated 
Gaussian random variables, and Gaussian random variables are uncorrelated if and only 
if they are independent. Thus, X(t) and X(t) are independent for the stationary Gaussian 
situation. The integrals can be evaluated explicitly, giving 


Note that the ratio g x I g x is exactly the co c characteristic frequency defined in 
connection with the energy-based amplitude A 2 in the preceding section. Also note that 
the maximum value of the crossing rate occurs when u = 0 and it is simply 
vj^(O) ~ co c /(2k). Thus, the rate of upcrossings of zero by a mean-zero stationary 
Gaussian process is simply the energy-based average frequency divided by 2k. The 
factor of 2k, of course, comes from the fact that co c represents a frequency in radians per 
second, while the rate of mean-upcrossings represents a frequency in cycles per second, 
or Hz. This relationship between mean-crossing rate and characteristic frequency could 
have been anticipated for any narrowband process, but it holds exactly for a stationary 
Gaussian process regardless of the bandwidth. 

Since a peak of X(t) occurs whenever X(t) = Oand X(t)< 0, the rate of 
occurrence of peaks of (X(f)} is exactly the rate of downcrossings of the level zero by 


Similarly, the rate of occurrence of valleys of {X(t)} is v v (t) = v' x ((),t) . For a stationary 
Gaussian process v P = v v = g x I (2kg x ) , just as v + x (0) = g x I(2kg x ) . 

Knowledge of the rates of occurrence of peaks and crossings can also be used to 
define a measure of bandwidth of a stochastic process. In particular, any sufficiently long 



2 KG x 


(31) 


v P ( t ) = v}(0 ,t) = J°^ | w | Px(t\x(t) (°> w ) dw = Px(t) (O)]^ I w I Px(t) H X (0 = 0 )dw (32) 
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continuous time history of a process must have at least as many peaks as it has 
upcrossings of any level, and the number of peaks is expected to be only slightly larger 
than the number of upcrossings of zero for a mean-zero narrowband process. Thus, the 
irregularity factor, defined as IF = v x (0,t)/v P (t) can be used as a bandwidth measure for 
any stochastic process. Its range of possible values is from zero to unity, and it tends to 
unity for a narrowband process, just as for the previously defined bandwidth parameters. 

For a stationary Gaussian process the irregularity factor becomes IF = /{<j x o x ) . 

Somewhat surprisingly, IF for this class of processes is exactly the same as the a 2 
bandwidth parameter defined in (26). This allows estimation of its value from occurrence 
rates for crossings and peaks of a time history, as well as from integrals in the frequency 
domain, provided that a % is finite and a 2 P 0. 


PROBABILITY DISTRIBUTION OF PEAKS 


The probability distribution of the peaks (local maxima) of {X(t)} can be found by 
a procedure that is basically the same as that used in deriving the rates of occurrence of 
crossings or peaks. An occurrence rate v P [t;X(t) < u] for peaks below a level u is defined 
as /’(peak < u during [tJ + Af\)IAt just as v P (t) = /’(peak during [t,t + At]) /At. What 
is commonly called the distribution of peaks is really the conditional distribution given 
the existence of a peak at that time: 


Fp(t )( u ) - 


/’(peak < u during [t.t + At]) _ v P [t; X ( t ) < u\ 
/’(peak during [t, t + At]) Vp(t) 


LL 1 2 1 Px{t)X(t)X(t)^^ dwdz 


(33) 


and taking a derivative with respect to u gives the probability density function as 


Pp(t)( u ) 


Ll Z I Px(t)X(t)X(t)( U ’®’ Z W z 

\_J z \Px(t)X(t)(°’ z)dz 


(34) 


Either (33) or (34) describes the probability distribution of any peak that occurs 
within the vicinity of time t. The probability that the peak is within any given interval can 
be found directly from (33) or from integration of (34), and these formulas are also 
convenient for evaluating other quantities such as the mean value and variance of P(t) . 

Note that the probability distribution of the peak P(t) depends on the joint 
probability distribution of X(t), X{t), and X(t ) . This joint distribution is relatively 
simple for the special case of a stationary mean-zero Gaussian process. All that is needed 
is the covariance matrix of the three random variables, and this involves only the three 
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standard deviations and the correlation coefficient between X(t) and X(t), since X(t) is 
uncorrelated with the pair [X(t),X(t)] for any stationary process. It is easily shown, 
however, that the correlation coefficient between X(t) and X(t) is exactly the negative 
of the a 2 bandwidth parameter. Thus, this bandwidth parameter plays an important role 
in the distribution of the peaks of {X(t )} . The resulting formulas are 


and 


^(M )= 2 g 

(2n ) <j x 


F P (u) = d 


/[2al-(l-«2 2 )] i c -u 2 /(2al) 


0 


<J 


(1 ~ai) 1/2 a x ) 


-ffig 




-u 2 /(2* 2 x ) 


a 2 u 


0 


V{\-al) m a x ) 


a 2 u 


U1 -al) m cj x ) 


(35) 


(36) 


1/2 

in which 0(y) = [1 + erf(v/2 cry)]/ 2 is the cumulative distribution function F Y (y ) for a 

mean-zero random variable Y. The formulas in (35) and (36) are commonly referred to as 
the Rice distribution, in recognition of Rice’s pioneering work on this problem [1945], 

The limiting forms of this distribution for a 2 = 1 and a 2 -0 yield important results 
regarding the peak distribution. For the narrowband situation with a 2 approaching unity, 
some of the arguments in pp(u) and Fp(u) tend to infinity and the limit for a 2 =1 is 
exactly the Rayleigh distribution of (27) that describes both the A x (t) and A 2 (t) 
amplitudes of the Gaussian process. The agreement of the peak distribution and the 
amplitude distribution of the limiting narrowband process is consistent with the fact that a 
narrowband process can be considered to be a hannonic function with slowly varying 
amplitude and phase. Because the narrowband amplitude varies slowly, each peak of the 
narrowband process is equal to the amplitude of the process at that instant of time, so it is 
not surprising that the two quantities have the same probability distribution. For the 
opposite extreme situation with a 2 = 0 , the probability distribution of the peaks becomes 
exactly the Gaussian distribution of X{ t ) , which is consistent with the fact that the rate of 
occurrence of peaks is infinite for a 2 - 0 . That is, there may be peaks everywhere along 
the process so the distribution of peaks is same as the distribution of X(t). 

A convenient feature of the Rice distribution in (35) and (36) is that it also describes 

9 1 /O 

a random variable R = a 2 R x +(l-a 2 ) ~R 2 , in which R { and R 2 are independent 
random variables with R j having the Rayleigh distribution of the A(t) amplitude and R 2 
having the mean-zero Gaussian distribution of X(t) [Krenk, 1978]. 

The distribution of peaks of a mean-zero narrowband process is sometimes 
approximated by a function that can be obtained from knowledge only of the crossing 
rates of X(t). The rationale is that in the narrowband case one can ignore the possibility 
of peaks below zero or valleys above zero. With this simplification, a peak occurs within 
the interval [u,u + Au] if and only if there is an upcrossing of the level u that is not 
followed by an upcrossing of the level u + Au. This approximation then implies that the 
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expected number of peaks within an interval [u,u + Au] is the difference between the 
number of upcrossings of the level u and the number of upcrossings of the level u + Au, 

and the expected rate of occurrence of peaks in the interval is v + x (u) - v x ( u + A a). After 
nonnalization this gives 


- 1 rfl/ t(") (37) 

v x(Mx) dll 

for u> 0. For the special case of a stationary Gaussian process this result is exactly the 
Rayleigh distribution, which is the true answer only in the limiting case with a 2 = 1. For 
a non-Gaussian process, as well, one can anticipate that (37) will be asymptotically 
correct for a 2 approaching unity. The approximation may be significantly inaccurate for 
some more broadband situations. 



FIRST-PASSAGE AND GLOBAL MAXIMUM 

Finding the probability distribution of a global maximum within an interval [0,t] is 
considerably more complicated than finding the distribution for a peak (i.e., a local 
maximum). The extreme value distribution for the process { X(t}} and any fixed time t is 
defined to be the distribution of the random variable 

X(t) = max X(s) (38) 

0 <s<t 


which makes {X(t)} a stochastic process. Note that even for a stationary { X(t)} , one 
must expect that { X(t)} will be nonstationary, since larger and larger values of X(t) will 
generally occur if the period of observation is extended. Letting L x {u,t ) denote the 
cumulative distribution function of X(t) gives 

L x(u,t) = F (w) = P\X(t) < u] sf[I(«) < u: 0 < s < t~\ (39) 

in which the notation on the final term means that the X(s) < u inequality holds for all 
the given 5 values. This L x {u,t) function is sometimes called the probability of survival, 
which is certainly appropriate if u denotes a critical value for \X(t)) corresponding to 
some failure mode of the system. The probability density function for the extreme value 
is simply the derivative dL x {u,t)l du, and from this information one can also calculate 
the mean and variance of the extreme value. 

An alternative problem that is almost equivalent to extreme value analysis involves 
the random quantity called first-passage time. Let T x (u ) denote the first time (after 
time zero) at which X(t) has an upcrossing of the level u. That is, 
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X[T x (u)] = u, X[T x (u)]> 0, and there has been no upcrossing in the interval 
0 < t < T x (u). This T x (u) quantity is a random variable for any given u value. 

It should be noted that some caution is necessary in analyzing either the global 
maximum response or the first-passage time for a process |7(t)} that has a nonstationary 
mean- value function ju Y (t). In this situation the event { Y(s) < v: 0 < .s' < /} is 
\X(s) < v + jl/ y ( s ) ; 0 ^ s < (} for the mean-zero process X(t) = Y(t)~ ju Y (t). Similarly a 
crossing of the level v by 7(7) coincides with a crossing of v-ju Y (t ) by X(t). This 
amounts to a problem with a variable barrier level, and such problems are not explicitly 
studied here since they are less frequently encountered in practice. 

To see the relationship between the first-passage time and the extreme value 
distribution consider the event \X(s) <u:0<s<t} that appears in (39). This event can 
also be written as {X(0 )<u,T x (u)>t}, since X(s) can be less than u throughout the 
time interval only if it starts below u and does not have an upcrossing during the time 
interval. Thus 

LxM = TOO) < u\P[T x (u) > 1 1 X(0) < u\ 

= L x (u,0)P[T x (u)>t\X(0)<u] 
and 

ptM<\x(o : < 41 > 

Thus, the partial derivative of L x (u,t ) with respect to t gives the conditional probability 
density of the first-passage time, just as the partial derivative with respect to u gives the 
probability density of the extreme response X(t). This close relationship between the 
extreme value problem and the first-passage problem is not always mentioned in the 
literature, with some authors using only one terminology and some using only the other. 

In some problems, one can neglect the conditioning in (40) and (41), treating T x {u) 
as independent of X(0) . In particular, this is true in problems with an initial condition of 
F\X({)) <u\ = 1 , such as when the system is known to start at X(0) = 0 . Even in situations 
in which there is no specific information of independence, it is usually true that the effect 
of the initial-value conditioning is significant only for small values of t. 

It is often convenient to write the probability of survival in an exponential form of 
L x (u, t ) =L x (u, 0) cxp[- | o ? )) x ( u,s ) dsj (42) 

which follows directly from a definition of r/ x (u,t) =\-\l L x {u,t)\dL x {u,i)ldt . This 
gives 

rj x (u,t)At =T , (upcrossing in [t,t + At] X(0) < n,no upcrossing prior to t ) (43) 
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Comparing (43) with the upcrossing rate in (28) shows that r/ x (u,t) is a conditional form 

of v x (u,t). Thus, r/ x (u,t) can be considered to be the conditional rate of upcrossings of 
the level u, given the initial condition and the fact that there has been no prior upcrossing. 
The idea of a first-crossing rate seems contradictory, though, since there can only be one 
first-crossing for a given time history. If exceeding the level u is considered to 
correspond to a failure of the system, then r/ x (u,t) is what is called the hazard function in 
reliability theory. 

Unfortunately, it is not easy to calculate r/ x (u,t) in order to find the probability 
distribution of either the first-passage time T x (u) or the extreme value X(t). The only 
rigorous, general relationship for r/ x (u,t ) is the so-called inclusion-exclusion series, 
which involves an infinite sequence of rather complicated integrals involving the joint 
distribution of X and X at multiple time values [e.g., see Madsen et ah, 1986]. 

In many practical situations r/ x (u,t) tends asymptotically to a function r/ x (u) that 
is independent of t as t becomes large. In particular, this asymptotic behavior is likely if 
|X(t)} is stationary, or becomes stationary for large values of t. The reason is that most 
physical processes have only a finite memory, in the sense that conditioning by an event 
far in the past has little effect on the current behavior of the process. When this is true 
one can use an approximation 

L x (u,t) « L 0 g - 'fr * for large t (44) 


If there exists a time value T 0 such that r/ x (u,t) « r/ x (u) for t>T 0 , then (44) 
agrees with (42) for t > T 0 if 

L 0 = L x (u,0)QxJ-\^jj x (u,s)-rj x (u)]ds] (45) 


The value of L {] then depends on the way in which r/ x (u,t) for small t differs from 
r/ x (u). For example, if (X(t)} is stationary for all time t, then rj x (u,t) for small time is 
usually greater than the stationary value r/ x (u) and L 0 is smaller than L x (u, 0)< 1. On 
the other hand, a common nonstationary process has the “zero-start” condition of 
X(0) = 0 and X(0) = 0 , as applies to the response of a dynamical system that starts from 
a condition of rest. In this case, rj x (u,t) starts from zero and increases as the response 
grows, giving L 0 > L x (u, 0) = 1 . These behaviors are illustrated in Figure 4. 

Note that (42) and (45) divide the probability of failure into two parts. The 
exponential tenns relate to the probability of crossing the level u, while the multiplier 
L x (u, 0) accounts for the fact that a process may start above the level u, which is given 
by P\X(Q) < u\. For the mean-zero narrowband process, note that there is a high 
probability of an early upcrossing if A( 0)>u. Since the narrowband amplitude varies 
slowly with time, it is very likely that any X(t) sample with 4(0) > u will exceed u 
during the first cycle. One way to approximate this effect is to treat the event 4(0) > u as 
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Figure 4. Effect of initial conditions on Lq in (44). 

though it were failure at time zero, so that the multiplier L x (u, 0) in (42) and (45) is taken 
as P[A(Q)<u\. Clearly, this approximation does not accurately describe the reliability 
during the first cycle of the mean-zero narrowband process, but it can give useful results 
for later times, which are generally of greater interest. 

Any error in estimating L 0 in (44) results in a proportional error in the predicted 

reliability L x (u,t). An error in ij x (u), though, is much more important since it leads to 
an error in L x (u,t ) that grows exponentially. Thus, there has been considerable effort 
invested is estimating this limiting ij x (u ) function. The following section summarizes 
some of the results obtained by direct consideration of the first-passage time. 


FIRST-PASSAGE FORMULATIONS 

The simplest approximation of ij x (u,t ) is related to the similarity between (43) and 
(28). Since r/ x (u,t) is a conditional fonn of the upcrossing rate, simply neglecting the 
conditioning gives ij x {u,t ) « v x (u,t). For any process this gives 

L x (u,t) «Z x (w,0)exp^- J ^(u^dsj (46) 

and in the particular case of a stationary mean-zero process it reduces to the form in (44) 
with L 0 =L x (u, 0) which may be taken as either P[X(ti) < u] or P[A(0) < u\. The 
formula in (46) is commonly called the Poisson approximation of the first-passage 
problem. This name comes from the fact that if the integer-valued process that counts the 
number of upcrossings by X(t) were a Poisson process, then the crossing rate would be 
independent of the past history of the process and the conditioning neglected in (46) 
would be irrelevant. 

For a stationary |X(t)} process, the Poisson approximation gives the first passage 

time T x (u) as having a simple exponential distribution with P[T x {u)<t\ = l-e~ Vx ^ t . 
This then gives the mean and standard deviation of T x (u ) as both being equal to 
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[v^{u)] 1 . In the particular case of a stationary mean-zero Gaussian process, the crossing 
rate in (3 1 ) then gives 


L x (u,t ) « L x (u, 0)exp 


f ~ 0) c t e -u 2 l(2a 2 x ) 

v 2k 


m 


in which co c is again the characteristic frequency defined in connection with the energy- 
based amplitude A 2 . Also, the value of P[X(0) < u] or P[A{ 0) < u] for use as L x (u, 0) is 
easily calculated for the Gaussian process, for which A has the Rayleigh distribution. It 
should be noted that (47) can only be considered to be an approximation of L x (u,t) for 

u> 0. In particular, P[X(t) < u] = L x (u,t) must be an increasing function of u, and this 
condition is not met by (47) for u< 0. Thus, one must consider (47) to give a finite 
probability that X(t) is equal to zero: P[X(t) = 0] = T x ( 0 , 0 )e - AM 2 ^)_ This discontinuity 
in F'.Au) decays toward zero as t increases. 

The error of the Poisson approximation is most serious when the (X(t)} process is 
very narrowband. In that situation, an upcrossing of level u at time t is very likely to be 
associated with another upcrossing approximately one period later, due to the slowly 
varying amplitude of (X(t)} . Such a “clumping” of the upcrossing times is inconsistent 
with the Poisson approximation that the times between upcrossings are independent. On 
the other hand, when u is very large it is found that the independence assumption seems 
to be better. There are few general results that apply to all (X(f)} processes, but for 
stationary Gaussian processes it has been demonstrated that rj x (u,t) does tend 
asymptotically to v x (u) as u tends to infinity [Cramer 1966]. Thus, the Poisson 
approximation is best when the {X(t)} process is very broadband and/or the level u is 
very large. In some narrowband situations, it may be significantly in error. 

From the fonn of (42) it is clear that an overestimation of r/ x (u,t) results in an 
underestimation of L x (u,t). An error of this type is usually considered to be 
conservative, because it overestimates the probability of failure due to large excursions. 
Furthermore, in most practical situations it is true that r/ x (u,t)< v x (u,t), so that the 
Poisson approximation does underestimate L x (u,t). Some caution is appropriate, though, 
since there are situations in which r/ x (u,t)> v x (u,t). In particular, if the level u is so 
small that P\X(t) < it] is also small and if there is an initial condition that A(0) < u, then 
it is very likely that X(t) will quickly have an upcrossing of u. Mathematically this 
requires that L x {u,t ) approach zero as it -oo, for any finite t value. This is ensured for 
any choice of L x {u, 0) only if r/ x (u,t ) tends to infinity for w— »- oo, at least for t~ 0, 
clearly requiring that rj x (u,t)> v x (u,t). This small-r/ situation, though, is not often of 
practical importance in reliability analysis. 

There are several ways in which the Poisson estimate has been modified to give 
somewhat better results. One of the simplest is an attempt to address, in an approximate 
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way, the noted underestimation of r/ x (u,t) for situations in which u is so low that the 
average time between upcrossings may be much larger than the time until the first 
upcrossing. Instead of the usual Poisson assumption that the probability distribution of 
the first-passage time T x {u) for a stationary process is the same as that of the total time 
interval between successive upcrossings, the modification is based on the idea that T x (u) 
has more in common with the time spent below the level u between successive 
upcrossings. This reduces the estimate of the time until first passage and gives more 
conservative estimates of failure probabilities for small values of u but has little effect for 
large u values. A simple approximation of this time increment below u gives 
r/ x (u,t ) « v x (uj)/ F X(l) (u) [Lutes and Sarkani, 2004]. This approximation does give the 
desired behavior that L x (u,t) tends to zero as u tends to negative infinity. This same 
approximation of ij x {u,t ) was earlier obtained by Ditlevsen [1986] by a somewhat 
different method of reasoning. 

Attention will now be focused on more significant modifications to the Poisson 
approximation. In particular, the aim of these approaches is to give better results for 
rj x (u) for a narrowband process, in which the upcrossing times cannot be independent 
because of the slowly varying amplitude { A(t )} . The simplest extreme value 
approximation based on the slowly varying amplitude process amounts to assuming that 
the extreme value of the mean-zero (X(t)} process is the same as that of {A(t)j , and that 
the Poisson approximation applies to {A{t )} . This gives L x (u,t) « L A (u,t), in which 
L A (u,t) is the same as (46) except that v x (u) is replaced by v A (it). By exactly the same 
reasoning as was used in suggesting that v x (u)/F x ^ t ^(u) was preferable to v x (u) for the 
original Poisson approximation, one can also say that v A (u)/F A ^(u) is preferable to 
v A (it) for this Poisson-amplitude approximation. 

For the special case of a stationary mean-zero Gaussian process, the v A (u) crossing 
rate has been evaluated for both the Cramer/Leadbetter and energy-based amplitude [e.g., 
Lutes and Sarkani, 2004]. The results are 
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and the probability estimate for first-cycle crossing of a narrowband process is 
F A ( ( )( u ) = 1 -e~ u /{2(Jx \ Thus, implementation of the Poisson-amplitude approximation 
is very straightforward. 
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It is generally presumed that the Poisson assumption is much better for {A{t)} than 
for {X(t )\ , and that the amplitude-crossing approach is conservative, since the maximum 
of X(t ) cannot exceed the maximum of A(t). Numerical studies show that the Poisson- 
amplitude approximation does give significantly improved estimates of r/ x (u,t) for small 
to moderate values of u values, but it is seriously in error when u is large. Whereas the 
original Poisson approximation of rj x ( u)»v x (u) becomes asymptotically correct as u 

goes to infinity, v^(u) / v x (u) grows without limit as u increases, giving grossly overly 
conservative estimates of ?l x (u). 

An improved approximation of the extreme value distribution can be obtained by 
estimating the fraction of the upcrossings by A(t) that are accompanied by upcrossings 
by X(t). The i] x (uj) conditional crossing rate can then be taken to represent the rate of 
occurrence of this subset of the amplitude-upcrossings. Vanmarcke introduced such a 
scheme [Vanmarcke, 1972; Corotis et ah, 1972; Vanmarcke, 1975] and Madsen et al. 
[1986] derived a very similar result. Some of the assumptions in these approximations 
will be summarized here, along with the fonnulas. 

Let the random variable 7] denote the time between an upcrossing of u by A(t) and 
the subsequent downcrossing by A(t). Then 7] represents the duration of an interval with 
A(t) > u. If T\ is large, then it seems almost certain that X(t) will have an upcrossing of 
u within the interval, but if 7] is small, then it seems quite likely that no upcrossing by 
X(t) will occur. Vanmarcke approximated this relationship by 


P[no upcrossing by X(t)\ 7] = r] « [1 - v + x (0,t) r]C/[l- v x (0,t) r] (50) 


Considering [v£(0,i)] _1 to represent the period of an average cycle of the (X(t)} 
process, this approximation amounts to saying that an upcrossing by X(t) is sure if 7] 
exceeds the period, and the probability of its occurrence grows linearly with 7] for 7] 
less than the period. Even though this approximation is crude, it is substantially better 
than simply assuming that an upcrossing by X(t) occurs in connection with every 
upcrossing by A(t). To calculate the unconditional probability of an upcrossing in the 7] 
interval, it is necessary to have a probability distribution for 7] , and this is assumed to be 
the exponential distribution p T ^ ( r) = X T ,:[T ' V7s(7]). Based on the usual formulation that 

the mean time between upcrossings is [v\(u,t)Y , the mean time below u during the 
interval is estimated as 7?(7]) « P[A(t) > u\/v^(u,t). Using this along with 
r/ x (u,t ) « if[(i/,t)f’[upcrossing by X{t) during 7]] gives 


r/ x {u,t) ~ P[A(t) > u] v x (0 ,t) 


1-exp 


\_P[A(t)>u]v + x {Q,t)\ 


(51) 
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In the limit as u tends to zero, P[A(t)>u\ tends to unity and >] x (uJ) tends to 
v\{u,t), as when considering each upcrossing by A(t) to correspond to a crossing by 
X(t). For large values of u, P[A(t)> u] is very small, so v\{u,t) » P[A(t) > u\ v x (0,t) 
and (5 1) gives rj x (u,t) « P[A(t) > it] n^(0,t) . For the special case of a Gaussian process 
it can be shown that this latter limit is identical to v x (u,t), so the approximation agrees 
with the results from the assumption of Poisson crossings by X(t). For a non-Gaussian 
process, these two results for large u values may not be identical, although they are 
expected to be quite similar. 


As in our other approximations, one can expect to obtain better results for small u 
values by including the effect of the initial condition. Because this estimate of r/ x (u,t) is 
a modified version of the amplitude-crossing rate, it is reasonable to do this by dividing 
by P[A(t)< u], giving 


rj x (u,t) 


P[A(t)> u]v x (0,t) 


P[A(t ) < u] 


1-exn 


-v + A (u,t) 


P[A(t)> u]v x (Q,t)_ 


(52) 


which agrees with v + 4 (u, t)/ F A(( ^(u) in the limit for u near zero. 


It should be noted that for a general non-Gaussian {A]T)} process, (52) is not 
identical to Vanmarcke’s result and that his derivation uses somewhat more sophisticated 
assumptions about the behavior of {X{T)} . For the special case of the Gaussian process, 

though, it can be shown that P[A x (t) > u\ v x (0,t) = v' x (u,t), so (52) with the Cramer and 
Leadbetter definition of amplitude does become identical to Vanmarcke’s form of 


T ! X (u,t) ~ v x {u,t) 


1 - exd 
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A t ( u A) 


v x (u,t) 
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(53) 


It can be expected that the two approximations will also give similar results for 
other processes that do not differ greatly from the Gaussian distribution. 

For the special case of a stationary, mean-zero, Gaussian |A(t)} process, 
Vanmarcke [1975] also offered an empirical correction that improves the approximation 
of the conditional crossing rate, at least for small to moderate values of u. For this 
situation, the ratio v^(u)/v x (u) of (53) is [2^(1 - a\)\ {l2 {ul <j x ) in which the term 
(1-aj 2 ) introduces the effect of the bandwidth of the process. As an empirical correction 
for effects not included in the derivation of these equations, Vanmarcke suggested 
replacing the (1 -a l ) term with (1 -a\) ' . This gives the modified Vanmarcke 
approximation as 
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In many problems, failure may occur due to large excursions of X(t) in either the 
positive or negative direction, whereas all the development up until now has been 
concerned only with the probability that X{t ) remains below +u. However, the event of 
X(t) remaining within \-u,+u\ is exactly the same as the event of \X(t) \ remaining 
below the level u, and this allows the new problem to be written in the form of (42), but 
with Z| X |(m, 0) and r/^(u,s) used in place of L x (u, 0) and r/ x (u,s), respectively. The 

tenns double-barrier problem and single-barrier problem are often used to distinguish 
between upcrossings by | X(t) | and X(t), respectively. Of course, one can also consider 
double-barrier problems in which the levels of interest are not symmetric. Modifications 
to the various first-passage approximations will be summarized for the symmetric 
double -barrier problem. 


The double-barrier Poisson approximation is that r/^(u,s) = v x (u,s)+ v x {-u,s). If 
the distribution of X{t ) and X(t) is symmetric, this then gives i] x (u,s) « 2 v x (u,s), 
indicating that the decay of L with increasing t is twice as fast as for the single-barrier 
problem. The small-time modification of dividing V\ X \(u,s) by F\ X \{u) has more effect 
than the corresponding adjustment for the single barrier problem, since F x (u) tends to 
zero as u tends to zero. The approximation of // by v' A (u,s) or v A (u,s)l F x (u) is equally 
valid for the double-barrier problem as for the single-barrier situation. That is, for the 
mean-zero process, the extreme distribution of A(t) can be considered an approximation 
of the extreme distribution of | X(t) |, just as well as for the extreme distribution of X(t). 
In fact, it seems likely that the approximation will be better for | X(t) | than for X(t). 


Adapting the Vanmarcke approximation to the double-barrier problem requires 
estimation of the probability that an upcrossing by A(t) is accompanied by an upcrossing 
by | X(t) j during the interval of length 7j between an upcrossing of u by A(t) and the 
subsequent downcrossing. This probability is estimated as zero for 7j = 0 and growing 
linearly to unity for 7j equal to the half-period [uj- (()./)] '/2. Other assumptions are the 
same as for the single-barrier situation and the result corresponding to (52) is 
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The corresponding modification of Vanmarcke’s formula in (53) is 
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Figure 5. Simulation and various approximations. 


which is identical to using the Cramer and Leadbetter amplitude in (55) for a mean-zero 
Gaussian {X(t)\ process. Using the empirical adjustment in (54) for a stationary process 
changes this to 
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An illustration is given in Figure 5 of the effects of the various assumptions on the 
double-barrier estimates of TJ\ X \(u) for a stationary, mean-zero, Gaussian, narrowband 

process with a x = 0.994. The results are presented as the ratio )], x (it)/ v' x (u ) , a 
nonnalization introduced by Crandall et al. [1966] and commonly used since then. Thus, 
the results for the assumption of Poisson crossings by X{t) appear as unity on this plot, 

even though they are inaccurate for such a narrowband process. The results for the v\{u) 
approximation are based on the Cramer/Leadbetter amplitude. There are two curves each 
for Poisson X crossings and Poisson A crossings, showing the increase in the value of 77 
when the crossing rates are divided by Fj X |(w) and F A ( u ) , respectively. Only the results 
for Vanmarcke’s approximation give all the desired tendencies of rj x {u ): tending to 
infinity for u going to zero, being less than v x (u) for intennediate values of u, and 
tending to v x (u) for large values of a. 

The plot in Figure 5 also includes some simulation results for one particular 
stationary, mean-zero, Gaussian process with the specified value of a x [Lutes and 

Sarkani, 2004]. It is seen that the simulated values of >] x (u) are smaller than the 
predictions from any of the approximations, except when u is less than about 1.2 a x . 
Also, the simulation data clearly show that the rj x (u)/ v x (u) ratio has a minimum value 
when u is approximately 2a x . The analytical form due to Vanmarcke comes closest to 
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fitting the simulation data. In particular, the modified Vamnarcke form of (57) gives a 
minimum value of r]\ X \(u)l v x (u) at approximately the right u value, and the values it 

gives for the ratio in this vicinity are better than those of any of the other approximations 
that have reasonable behavior for smaller u values. Nonetheless, it must be noted that 
there is sometimes a significant discrepancy between Vamnarcke ’s formula and the 
simulation data. For u = 2cr x , for example, the approximation is about 70% above the 
value of 0.10 or 0.1 1 obtained from simulation, even though the modified approximation 
agrees almost perfectly with the data point from simulation for u = 4 a x • Any 
overestimation of r/^(u) gives L x (u,t) values that decay more rapidly with increasing t 

than do the values from simulation. For u = 2a x and large values of time, Vanmarcke’s 
formula will significantly overpredict the probability that |X(t)| has ever reached the 
level u, but this discrepancy is much smaller than that of the commonly used assumption 
of Poisson crossings by X(t). Also, the discrepancy is smaller when u is large, and this is 
often the region of primary interest in predicting failure. It should also be noted that the 
simulation data suggest the possibility that the modified Vamnarcke approximation in 
(57) may be nonconservative for u values greater than 4<j x . That is, the empirical 
adjustment in (54) consistently improves the results for moderate values of u, but it may 
be excessive for very large values of u. 


GLOBAL EXTREME FORMULATIONS 

Recall that X{t) was defined in (38) as the maximum value of the process (V(t)} 
within the interval [0,t]. There is an obvious relationship between the problem of 
describing this global maximum and the classical problem of describing the maximum of 
a set of random variables: R = max{R 1 ,R 2 ,---,R, ; } • The simplest version of this problem 
is when the Rj random variables are independent and identically distributed, which gives 

the distribution of R as 

F r (m) = Fr (m) P r (u) = n F R ~ l ( u)p R ( u ) (58) 

The fact that F R (u) must tend from zero to unity as u increases, means that 
F r ~\u) has a similar fonn, but is shifted further and further to the right along the u axis 

as n increases. Thus, the probability distribution of R for very large n values depends 
only on the behavior of the right-hand tail of p R (u). 

One way to apply this random variable idea to the stochastic process problem is to 
note that X(t) is the maximum peak of { X(t)\ during the interval [(),/] . For the mean- 
zero Gaussian process, ® c /(2;ra 2 ) is the expected peak-occurrence rate, with 
co c =<j x /<j x . Thus, the expected number of peaks is n P = (a> c t)/(2 7ra 2 ). If one now 
assumes that these n P peaks are all independent random variables then 

L x (u,t) = F- X ^^u) « F p p ( u ) = exp(<z) c tln[Fp(i/)]/[2^a 2 ]) (59) 
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which is exactly the same as the asymptotic crossing-rate fonn in (44) with Z 0 = 1 and 
r/x(u) = [-® c /( 2 ^ 2 )] HFp(u)] . 

It has been shown that the large-// results of the Poisson-crossings approximation 
can also be obtained by assuming that all the peaks within the [0,/] interval are 
independent random variables [Crandall, 1970]. To see this, recall that the Rice 
distribution of the peak P of a Gaussian process, as given in (36) is exactly the same as 
the weighted sum of two independent random variables, one having the Rayleigh 
distribution of the amplitude of { X(t)) and the other having the Gaussian distribution of 
{ X{t )} . The weighting factors for the Rayleigh and Gaussian components, respectively, 
are the bandwidth parameter a 2 and (1 -a 2 ) . The probability density of the Rayleigh 

distribution, as given in (27), is proportional to ue~ u /<2ax - 1 whereas that of the Gaussian 

distribution is proportional to e~ u /<2(7x \ From this comparison it is clear that the Rice 
distribution will be dominated by the Rayleigh component when u is very large. Since 
large-// is the situation of interest when n is large, this gives F P (u) from (36) tending to 

1 -a 2 e~ u 1(2(7 x ( and \n[F P (u)\ tending to {-a 2 e~ u /( 2(7 x ( ) . The asymptotic first- 

crossing rate is then r/ x (u) = (co c /27i)e~ u 1(2(7 x ) = v x (u) , which is exactly the same as 
the result of the Poisson-crossings assumption for X(t). 

For a narrowband Gaussian process there is an even simpler approach that gives 
asymptotic equivalence between (58) and the Poisson X-crossing approximation. 
Consider the interval [0,/] to be divided into n c cycles of period (2 7zlco c ), giving 
n c =(co c t)/( In'). Assume that the amplitude A(t) can be considered to be constant over 
each of these cycles and independent in different cycles so that (59) applies, giving 

L x (u,t) « F^ c (//) = exp(/u c tln[iyj (w)]/[2^]) (60) 

For the Rayleigh random variables, ln[F 4 ^)(//)] = ln[l-e _ “ 2(2(7 x > ] and for large // this 

tends to ( -e~ u j{2(7x ) ) , Again the asymptotic behavior is the same as the Poisson- 
crossings result with r/ x (u) = v x (u) . The formulas in (59) and (60) both give 
rj x (u ) > v x (i<) for small values of u. In fact, the curves are generally similar to those for 
?Ix(u)=v x (u)IF x (u). 

Naess and Gaidai [2008] recently introduced another approximation that may be 
viewed as a modification of the formula in (58). The approximation is based on the peaks 
of the time history, as in (59), but including conditioning based on the values of prior 
peaks. In particular, the approximation is that 

L x (u,t)«F^- k \u\B k ) (61) 

1 n 
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in which the conditioning event for peak P n is B k = {P n _\ < u,P n _ 2 < u,---,P n _ k < u) (i.e., 
the preceding k peaks have all been below u). The method used by Naess and Gaidai to 
estimate the critical quantity F P (u \ B k ) involves simulation of one or more relatively 

long time histories of X(t). It is easy to create a computer program to count peaks in a 
simulated time history and detennine the fraction of peaks above u for which there are k 
preceding peaks that are all below u. This fraction is used as the estimate of 
[1 - F P (u | B k )\ . One difficulty of such a simulation approach is that it may require a very 

long time history to find any significant number of peaks above a high u level. Naess and 
Gaidai, though, found that [1 -F P (u\B k )\ for large u is approximated by the form 

q 0 exp[-a(u-b) c ]. After evaluating the parameters to fit observed data, this form could 
be used to extrapolate [1 - F P „ (u \ B k )] to larger a values. The method has been 

demonstrated to give very good approximations of empirical data for k values of 3 or 4. It 
also has been demonstrated for a process with bimodal spectral density comprised of two 
narrow-band components. In this situation, though, the value of k must be large enough to 
include peaks over a selected number of lower-frequency periods. This can be a 
significant number when the frequencies of the components are well separated. 


A modified version of the largest cycle-amplitude approach is based on including 
some effect of correlation of the cycle amplitudes [Lutes, 2008]. The major simplification 
is to include only the effect of knowing that A(t- T c )< u when estimating the probability 
of A(t ) < u, in which z c is a representative period of the process. Obviously, this is very 
similar to the peak conditioning in (61) with k = 1 , but with the advantage of working 
with the distribution of amplitudes rather than peaks. The only new information needed is 


q(u) = P[A(t) < u | A(t 


* c) <“] = 777^ J*o J*o PA{t-r c )A{t)(y\Xl)- dv x dv 7 


F A ( U ) 


(62) 


or, in a form that is more accurate for numerical evaluations, 

?( M ) = p A{t-r c )A{t)(y i,v 2 )^iJv 2 ] (63) 

The result of this approximation is 

L x(u,t) = F ^ (ncTc) (u) * F A {u)q (ric ~ l) {u) = [F a ( u) Iq( u)] exp(7 In [q( u )] / r c ) (64) 


which gives the asymptotic first-crossing rate as ^ x (w) = -ln [q(u)\/ x c . Note that q(u) is 
bounded as F A ( u) < q(u) < 1, and it approaches the upper limit for a narrowband process. 
This gives an estimate of //^(w) that can be significantly smaller than that of the Poisson 
assumption for a narrowband process. 

Under certain circumstances the joint probability density for the amplitudes in (62) 
and (63) reduces to a relatively simple fonn. In particular, if the {A(0} process is 
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Gaussian so that the amplitude has the Rayleigh distribution, and if Gzz( t c)~Gxx( t c )’ 
for the {Z(t)} used in defining the amplitude, then the jointly Rayleigh distribution is 


PA(t-r)A(t)( u b u l) 


u \ u 2 

(1 ~P)Ox 


exp 


f 9 9 'N 

~( u \ + u i) 

\ 2(1 ~ P)c?xj 


p V ~U\U 2 
1(1 ~p)°x) 


(65) 


in which 7 0 (-) denotes the modified Bessel function of order zero. The new parameter p 
in (65) is the correlation coefficient of A (t) and A~(t+ n la> c ), so is always in the range 

9 9 a 

[0,1]. It can also be written as p = [G xx (z c ) + G^^x^l <j x . Although no closed-fonn 
solution has been found for the double integration in (62) or (63) for this probability 
distribution, numerical evaluation is relatively simple. For the Cramer/Leadbetter 
amplitude, G zz (z) = G xx ( z) for all values of z so that the Gzz( r c ) = G xx ( z c ) condition 
on (64) is automatically satisfied and one can use the natural choice of z c =2 n! <x> c . For 
the energy-based amplitude, the condition can sometimes be met by making a small 
adjustment in x c [Roberts, 1976]. 

For the Cramer/Leadbetter amplitude, the quantity G X7 { z) is the negative of the 
Hilbert transform of G xx ( z) and for the energy-based amplitude it is ~G' XX ( z) I co c . 
These quantities can be further simplified for the situation of a narrowband process. In 
particular, it is known that G^^z) is nearly hannonic with period z c for a narrowband 
process, and this nearly harmonic nature of the functions gives G zz (z c ) ~ G xx ( z c ) and 

G XZ ( r c) « Gjociz^. Thus, the distribution in (64) with p « G xx ( z c )/a x is a good 
approximation for any narrowband process. 

This method of considering conditioning of cycle-amplitudes as opposed to peaks, 
as in (64), has the advantage of avoiding numerical simulation of time histories. On the 
other hand, the method does require numerical integration of (63). In principle, the 
method can also be extended to include some amplitude-correlation over multiples of z c , 
giving something even closer to (61), but this would be significantly more difficult. 

The modifications needed in order to apply the methods in this section to the 
double -barrier problem are quite straightforward. The independent peak approximation 
of (59), now counts both peaks and valleys and half-cycles replace cycles in the cycle- 
amplitude approaches. This requires replacing n P with 2 n P in (59) and replacing n c with 
2 n c for cycle-amplitudes. The only change in the Naess and Gaidai approach is in the 
simulation to find the conditional distribution function in (61). 

Numerical results for the independent-amplitude method of (60) and the correlated- 
amplitude method of (64) are shown in Figure 6, along with the same simulation data as 
in Figure 5. The independent-peaks method of (59) and the correlated-peaks method of 
(61) are not useful for the process producing this simulation data, since the data are for a 
process with a 2 = 0 due to a x = cc . This limiting condition gives the peak distribution 

as Gaussian, so that it does not have the necessary tendency to the Rayleigh fonn for 
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Figure 6. Simulation and amplitude approximations. 


large u, as in all other situations. Figures 5 and 6 show that for this narrowband process 
the correlated-amplitude approach gives // values that are quite similar to those of the 
Vanmarcke method, but are somewhat more conservative for large u values. 


NUMERICAL RESULTS 

The significance of the preceding analytical results will now be illustrated by 
presenting figures and tables demonstrating the probability of | X(s) | exceeding a level u 
at least once during the interval 0 < s < t : 

P e (u,t ) = 1 -L\ X \(u,t) = P[X(t ) > »] = P[Px( u ) < (66) 

in which, X(t) is the largest value of | X(s) | within the set 0 < s < t and T x (u) is the 
time at which |X(t)| first exceeds the level u.. These results will be limited to the 
situation in which it is given that the process begins with | X(0) |< a and initial conditions 
will be ignored, so that (44) gives L\ X \(u,t) = exp[-rj\ X \(u)t]. 

The most comprehensive results will be presented for the Poisson approximation of 
the double-barrier first-passage probability, since it is relatively simple and is 
asymptotically correct for large u values. In particular, letting r/\ X \(u) = 2v x (u) and 
solving for time gives 


^^Xn[\-P e (u,t)]P K2ah (67) 

In 2 

Note that the left-hand side of (67) is a normalized time that can be taken as the number 
of “cycles” (or upcrossings of zero) of the mean-zero Gaussian process during the time 
interval [0,t]. Figure 7 then shows this cycle-count versus P e for several integer values of 


30 




u/a x . This can be viewed as the time required for P[X(t)> u ] to reach level P e for a 
stationary process or, equivalently, as the duration of a finite time history for which 
P[X(t) > u] = P e . As would be expected, the time value grows with increasing u and also 
with increasing P e . What is more interesting is the magnitude of the numerical values. 
For the 5-sigma level, for example, there is a 0.1% probability of exceedance within a 
time history of about 134 cycles, increasing to 1% probability for 1,350 cycles, and to 
50% probability for 93,000 cycles. Thus, there may be many situations in which 
exceeding 5 -sigma should not be considered to be a very rare event. Similarly, the 
exceedance of 3 -sigma may be considered to be quite common, having a 50% probability 
of occurrence within a time history of only 3 1 cycles. 


The likelihood of seemingly large extreme values can be illustrated in another way 
by using (67) to find the probability distribution of the largest extreme for a given value 
of t. In particular, the probability density of X(t) is simply the derivative of L x >(u,t) 
with respect to u, but with the modification that the Poisson approximation requires a 
finite value for L\ X \(fd,t), as discussed following (47). For the Poisson approximation of 
the double-barrier problem, this gives the probability density as 


p hl) (u) = e-°‘ ,l ’’m + 


r <»A 


V n J 




-exp 


-u 


\2<j x 


MS c ~u 2 /(2al) 

n ) 


( 68 ) 


This probability distribution can then be used to evaluate numerically the mean ( u v ) and 
the standard deviation ( <j - ) of the largest extreme within a time history of a given length. 

For illustration consider a time history of 1,350 cycles, for which it was previously found 
that there is a 1% probability of finding an extreme exceeding 5 a x . The probability 

density of X{t) for this situation is shown in Figure 8, and it is found numerically that 
jUj. « 4.03 <j x and cr^«0.31cr x . Thus, //^ + 3cr^ « 5<j x . For this particular example 


31 





there is a 1% probability of the largest extreme exceeding // - +3 cr~ — the mean plus 

three standard deviations of the random extreme. Again, it seems that 5<j x may be 
viewed as a value that is relatively likely to be exceeded in many studies of extreme 
values. 

Recall that the Vanmarcke approximations and the correlated-amplitude 
approximation all give results with L\ X \(u,t) decaying exponentially as ?] x <(u)t . Thus, 
any of these methods gives results similar to Figure 7 except that the time for P[X(t) > u] 
to reach level P e is increased by the ratio v x (u)h} x (u) . For the Vanmarcke 
approximations this ratio depends on the bandwidth parameter a 1 of X(t), and for the 

9 9 A 

correlated-amplitude approximation it depends on p = \G xx (n I co c ) + 

Tables 1 and 2 show the extent to which the time until exceedance is increased by 
use of the Vanmarcke approximations instead of the Poisson approximation. In particular, 
the numbers in the tables are multipliers to apply to the Poisson time values, as in (67) 
and Figure 7, to obtain a Vanmarcke estimate of the time for P\X(t)> it] to reach level 
P e for several values of u, and for several different narrowband processes. Note that 
Table 1 is for the original version of the Vanmarcke approximation, without the empirical 
correction factor introduced in (54) and (57), with this latter form being shown in Table 
2. Table 3 presents the same information for the approximation of correlated amplitudes, 
as in (64) and (65). 

The values of the parameters a 1 and p in the tables, of course, are related to the 
bandwidth of the process. To illustrate their significance, it may be helpful to consider a 
relatively well-known problem that does yield such narrowband processes as those 
presented here. In particular, consider the response of a single degree-of-freedom 
oscillator subjected to a white noise excitation, for which the bandwidth is governed by 
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the fraction of critical damping C, in the oscillator. For this example, the parameter a 1 is 
found to be 


<2i = 


(i-f 2 ) 1 ' 2 


2 2 J 

1 tan 1 


V 


n 




\\ 




(69) 


Similarly, the correlated-amplitude approximation is easily applied to this situation 
by using the A 2 definition of amplitude and choosing the time increment z c for one half- 

9 1/9 

cycle to be given by co 0 r c = n !{\ -£ ) , in which a> 0 is the undamped frequency of the 

oscillator. This gives 


p = exp[-2<^;r /(I - <^ 2 ) 1/2 ] (70) 

The results in (69) and (70) are also approximately correct for any broadband 
excitation of the simple oscillator, and they tend to a x ~\-2<^l n and p~\-2^n for 
very small values of 


Table 1 . Increase in time for Vanmarcke correction of (55) and (56) 




u = 3 

u = 4 

u = 5 

u = 6 

U = 1 

u = 8 

a 1 = 0.999 


6.39 

4.98 

4.09 

3.50 

3.08 

2.77 

a 1 = 0.99 


2.40 

1.97 

1.70 

1.53 

1.41 

1.32 

a 1 = 0.95 


1.43 

1.26 

1.16 

1.11 

1.07 

1.05 

a 1 = 0.9 


1.23 

1.13 

1.07 

1.04 

1.02 

1.01 


Table 2. Increase in time for Vanmarcke correction of (57) 




u = 3 

u = 4 

u = 5 

u = 6 

u = 7 

u = 8 

a 1 = 0.999 


11.45 

8.81 

7.16 

6.05 

5.26 

4.67 

a 1 = 0.99 


3.28 

2.63 

2.22 

1.95 

1.76 

1.62 

a 1 = 0.95 


1.63 

1.41 

Ml 

1.18 

1.13 

1.09 

a 1 = 0.9 


1.32 

1.19 

1.11 

1.07 

1.04 

1.03 


Table 3. Increase in time for correlated-amplitude correction of (64) 




u = 3 

u = 4 

u = 5 

u = 6 

u = 7 

u = 8 

p= 0.99 


8.27 

6.29 

5.05 

4.23 

3.65 

3.17 

35 

ii 

o 

M2 


2.64 

2.07 

1.72 

1.49 

1.34 

1.24 

p= 0.6 


1.38 

1.18 

1.08 

1.03 

1.01 

1.00 

p= 0.3 


1.11 

1.03 

1.01 

1.00 

1.00 

1.00 
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The values of a 1 in Tables 1 and 2 and the values of p in Table 3 have been chosen 
to correspond to similar values of £ in the single degree-of-freedom oscillator. In 
particular, the a x values of 0.999, 0.99, 0.95, and 0.9 correspond to C, values of 0.0016, 
0.016, 0.084, and 0.18, respectively, while the p values of 0.99, 0.9, 0.6, and 0.3 
correspond to C, values of 0.0016, 0.017, 0.081, and 0.19. The numbers in the first two 
rows of each table can be considered to be for very narrowband processes. Those in the 
first row ( £ = 0.00 1 6 ), in particular, are for a situation that is so narrowband as to be 
extremely rare in practice. 

The modified Vanmarcke results in Table 2, in particular, show that it is possible 
for this correction to give exceedance times that are an order of magnitude larger than 
those predicted by the Poisson approximation. The increases in exceedance time for the 
correlated-amplitude approximation in Table 3 are also significant for the same 
situations. However, these large increases only occur for very narrowband processes and 
a relatively low crossing level. Thus, the simple Poisson approximation is quite adequate 
for many practical situations. 


TRANSLATION PROCESSES 

There is one class of non-Gaussian processes for which it is particularly simple to 
apply all the results developed for the Gaussian situation. In particular, consider a non- 
Gaussian process {X(t)} that can be written as 

X{t) = g[Y{t)] (71) 

in which { Y(t)} is a normalized stationary Gaussian process and g(-) is a monotonically 
increasing function. The cumulative distribution of such a so-called translation process 
(A(t)} [Grigoriu, 1995] is then given by 

F x(t) [&(T)] = P{X(t) < g(y)] = P[Y(t)<y] = &{y) (72) 

Furthennore, the peaks of X{t) occur at the same time as those of Y(t) and the 
relationship between the peak magnitudes is Px = g(P Y ). In addition, the global maxima 
within a time interval [0,f] satisfy X(t) = g[F(/)], with a similar relationship for the 
global minima. Note that the global maximum and the global minimum of X(t) may 
have quite different probability distributions if g(- ) is not an odd function. 

In practice it is often useful to apply the translation procedure of (71) in order to 
approximate a non-Gaussian process {X(t)} by a translation process |g[T(0]} . One can 
always find an appropriate translation process by choosing g(-) to match the cumulative 
distribution function of g(Y(r)] to that of X(t). In particular, the inverses of the two 
distribution functions give u = F x ] [P(X < it)] and y = 0 1 [P(Y < v)] . Plotting u versus y 
for a range of probability values then gives the desired g(-) function. It should be noted, 
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though, that this procedure does not necessarily give {g J [X(f)]} as a Gaussian process. 
Rather, it assures that (g _1 [X(t)]} can be viewed as a collection of Gaussian random 
variables, but these random variables may or may not be jointly Gaussian, as is required 
for a Gaussian process. 

It is known that the peaks of the nonnalized Gaussian process (7(7)} are given by 
the Rice distribution 


P{Y<y) = d\ 


y 


(i -« 2 ) 


2x1/2 


\-a~,e 


- y 1 12 




«2 y 


j 


Ui-« 2 2 ) 1/2 ^ 


(73) 


The fact that g(v) is a monotonically increasing function assures that every peak of Y(t) 
maps into a peak of g[7(t)]. If {g \X(t)\ \ truly is a normalized Gaussian process, then 
the probability distribution of the peaks of g[7(7)] agrees with that of the peaks of X(t). 
Otherwise, {X(7)} is not truly a translation process. Of course, the Rice distribution of 
peaks is a necessary but not sufficient condition for Gaussianity, but it is one of the easier 
tests to apply in practice. 

Even in a situation in which the translation concept gives only an approximate 
description of the data being analyzed, it is fairly likely that the approximation will be of 
some value in predicting the probability of first-passage failure for a non-Gaussian 
process. 


SUMMARY 

The reliability problem can be formulated either as finding an estimate of the first- 
passage time T x {u) for a given value u of the stochastic process { X{t )} , or as estimating 

the largest value X(t) of {X(s) : 0 <s<t} for a given value of t. The probability 
distributions of these two random quantities are governed by the one function 
L x (u,t) = P[X(t)<u], since the probability density of X(t) is the //-derivative of 
L x {u,t ) and the probability density of T(a) is proportional to its /-derivative. The role of 
time is critical in the evaluations. For an ergodic unbounded process, such as a Gaussian 
process, the mean value of X{t) grows without limit as t becomes very large. 

One can always write L x (u,t) as an exponential function containing a time-integral 
of a quantity r/ x (t,u). For a stationary process, this new quantity depends on t only when 
t is small. Thus, it is efficient to concentrate on estimating the stationary value rj x {u ). 
For other than very small values of t this gives L x (u,t ) « Since L x {u,t ) can 

be viewed as a probability of survival, a conservative procedure is one that under- 
estimates L x (u,t), and for large u this requires an over-estimation of r/ x (u). On the same 
basis, it is desirable to under-estimate the value of L 0 . For a zero-start process that tends 
to a stationary condition this conservative estimate can be taken as L 0 «1. For a truly 
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stationary process a reasonable estimate is L 0 « P\A(t) < u] based on the amplitude of the 
process. This is more conservative than using P[X(t) < u] to estimate L 0 . 

The most commonly used method for approximating r/ x (u ) involves Poisson 
crossing times for X(t), and gives r/ x (u) as the same as the expected up-crossing rate 
v x (u). This approximation is conservative for all but small values of u, and is 
asymptotically correct for a Gaussian process as u becomes very large. A more 
conservative and better estimate for small u is obtained by using r/ x (u) « v x (u)/F x (u). 
The major shortcoming of the Poisson X-crossings approach is that its conservatism may 
be regarded as excessive for some narrowband processes for some values of u, with the 
greatest discrepancy occurring for u&2<j x . Using an alternate assumption of Poisson 

crossing times by A(t) leads to r/ x (u) « v + 4 (u) I F A (u) . This approximation provides a 
significant improvement for some values of u, but it can give the nonconservative result 
of r/ x (u) > v x (u) for large u values. A conservative estimate of r/ x (u) can be written as 
mm[v\(u) I F a (u), v x (u)], but this slightly complicates the computation of such 
quantities as the mean or standard deviation of X(t), since they involve integrations of 
the probability density of X{t) overall positive values of u. 

The Vanmarcke approximation of (52)-(57) provides the best available estimates of 
Tj x {u) using only analytical expressions. The formulas are significantly more 
complicated than for the Poisson crossing approaches, and they involve a factor a 1 that 
relates to the bandwidth of the spectral density of { X{t)} . Nonetheless, the method is 
easily implemented in a computer program. Similar, but slightly more conservative 
results are obtained from the correlated-amplitude approach of (62)-(65), which involves 
some straightforward numerical integration. 

The recently introduced correlated-peaks method of Naess and Gaidai seems to be 
capable of providing very good estimates of r/ x (u) for a variety of processes, including 
ones that are narrowband. Significant simulation of time histories of { X(t)} is required, 
though, to determine the values of key parameters in the formulation. 

Other analytical techniques summarized herein show the similarities and differences 
between the results of various assumptions, but none show improvements over the 
Poisson, Vanmarcke, and correlated-amplitude approximations. A plot and tables have 
been included to provide a simple overview of the length of time history for which 
P[X(t) > u] = P e for chosen values of u and P e for each of these approximations. 

The simple concept of a translation process has been introduced to provide an 
approximate method for applying the Gaussian techniques of first-passage analysis in a 
situation in which the process of interest is found to be non-Gaussian. 
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